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Abstract 

State-of-the-art ab initio techniques have been applied to compute the potential energy sur- 
face for the lithium atom interacting with the lithium hydride molecule in the Born-Oppenheimer 
approximation. The interaction potential was obtained using a combination of the explicitly cor- 
related unrestricted coupled-cluster method with single, double, and noniterative triple excitations 
[UCCSD(T)-F12] for the core-core and core-valence correlation and full configuration interaction 
for the valence-valence correlation. The potential energy surface has a global minimum 8743 cm -1 
deep if the Li-H bond length is held fixed at the monomer equilibrium distance or 8825 cm -1 
deep if it is allowed to vary. In order to evaluate the performance of the conventional CCSD(T) 
approach, calculations were carried out using correlation-consistent polarized valence X-tuple-zeta 
basis sets, with X ranging from 2 to 5, and a very large set of bond functions. Using simple two- 
point extrapolations based on the single-power laws X~ 2 and X -3 for the orbital basis sets, we 
were able to reproduce the CCSD(T)-F12 results for the characteristic points of the potential with 
an error of 0.49% at worst. The contribution beyond the CCSD(T)-F12 model, obtained from full 
configuration interaction (FCI) calculations for the valence-valence correlation, was shown to be 
very small, and the error bars on the potential were estimated. At linear LiH-Li geometries the 
ground-state potential shows an avoided crossing with an ion-pair potential. The energy difference 
between the ground-state and excited-state potentials at the avoided crossing is only 94 cm -1 . 
Using both adiabatic and diabatic pictures we analyse the interaction between the two potential 
energy surfaces and its possible impact on the collisional dynamics. When the LiH bond is allowed 
to vary, a seam of conical intersections appears at C*2v geometries. At the linear LiH-Li geometry, 
the conical intersection is at a Li-H distance which is only slightly larger than the monomer equilib- 
rium distance, but for nonlinear geometries it quickly shifts to Li-H distances that are well outside 
the classical turning points of the ground-state potential of LiH. This suggests that the conical 
intersection will have little impact on the dynamics of Li-LiH collisions at ultralow temperatures. 
Finally, the reaction channels for the exchange and insertion reactions are also analyzed, and found 
to be unimportant for the dynamics. 
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I. INTRODUCTION 



Ultracold molecules offer new opportunities for scientific exploration, including studies 
of molecular Bose-Einstein condensates, novel quantum phases, and ultracold chemistry. 
For molecular interactions that take place at microKelvin temperatures, even the smallest 
activation energy exceeds the available thermal energy. This opens up new possibilities for 
controlling the pathways of chemical reactions (see, e.g., Ref. 

A major objective of current experiments on cold molecules is to achieve quantum de- 
generacy, particularly for polar molecules. Two approaches are being pursued: indirect 
methods, in which molecules are formed from pre-cooled atomic gases, and direct methods, 
in which molecules are cooled from room temperature. There have been very substantial 
recent advances, particularly in indirect methods. In particular, the JILA Q and Innsbruck 
j^l groups have formed deeply bound ground-state molecules at temperatures below 1 /iK, 
by magnetoassociation of pairs of ultracold atoms followed by coherent state transfer with 
lasers. Methods that form ultracold molecules from ultracold atoms are however restricted 
at present to species formed from heavy alkali-metal atoms. 

n n 

Direct methods, such as buffer-gas cooling |4|, Stark deceleration |5j], crossed-beam col- 
lisional cooling [6] and Maxwell extraction 7j, are applicable to a much larger variety of 
chemically interesting molecules. However, these methods cannot yet reach temperatures 
below 10 to 100 niK. Finding a way to cool these molecules further, below 1 mK, is one of the 
biggest challenges facing the field. The most promising possibility is so-called sympathetic 
cooling, in which cold molecules are introduced into an ultracold atomic gas and thermalize 
with it. Sympathetic cooling has been successfully used to achieve Fermi degenerac y in 



6 Li [8fl and Bose-Einstein condensation in 41 K |9j, and for producing ultracold ions 10Nl2|. 
However, it has not yet been achieved for molecular systems, although there are theoretical 
proposals for experiments in which ultracold NH or ND3 molecules are obtained by collisions 



with a bath of colder atoms such as Rb, Mg or N 



13- 



The group at Imperial College London recently succeeded in producing samples of cold 

17] using Stark deceleration. LiH is 



LiH molecules in the first rotationally excited state [16 
an attractive molecule for cooling, since it has large dipole moment and light mass, so that it 
can be controlled easily with fields. It has a relatively large rotational constant (7.5 cm^ 1 ), 
which opens up the possibility of producing cold molecules in a single excited rotational state. 
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There is proposal to produce ultracold LiH molecules by sympathetic cooling with Li 18]. 
However, sympathetic cooling can be successful only if the rate of elastic (thermalization) 
collisions is large compared to the rate of inelastic (deexcitation) collisions, which cause trap 
loss. The main objects of the present paper are to explore the interaction between Li atoms 
and LiH molecules, to understand the nature of the interaction between these two species, 
and to obtain a detailed and accurate potential energy surface for the Li-LiH system. 

The results of scattering calculations at ultralow temperature are very sensitive to the 
details of the interaction potential 13|, [l5|. For systems containing heavy atoms, the meth- 



ods of quantum chemistry currently available cannot generate interaction potentials with 
accuracy better than a few percent. This limitation is caused by approximate treatments of 
correlation effects and relativistic contributions. With potential energy surfaces of moder- 
ate precision, it is usually possible to extract only qualitative information from low-energy 
collision calculations. By contrast, Li-LiH is a light system containing only 7 electrons 
and state-of-the-art ab initio electronic structure calculations can be performed with no 
significant approximations. It therefore offers a unique possibility to produce a very pre- 
cise interaction potential, which will allow a quantitative description of Li-LiH collision 
dynamics, even in the ultralow temperature regime. 

In electronic structure calculations one aims at approaching the exact solution of the 
Schrodinger equation, as closely as possible within the algebraic approximation. In practice, 
this is accomplished by combining hierarchies of one-electron and iV-electron expansions. 
The accuracy increases across the hierarchies in a systematic manner, allowing the errors in 
the calculations to be controlled and a systematic approach to the the exact solution to be 
achieved. The standard iV-electron hierarchy employed in electronic structure calculations 
consists of the Hartree-Fock (HF), second-order M0ller-Plesset perturbation theory (MP2), 
coupled-cluster with single and double excitations (CCSD), and coupled-cluster with single, 
double, and approximate noniterative triple excitations [CCSD(T)] models, with the latter 
recovering most of the correlation energy. Thus, CCSD(T) constitutes a robust and accu- 
rate computational tool nowadays. All these models are size-consistent, which means that 
the interaction potential shows the correct dissociation behaviour at large intermolecular 
distances. In contrast, methods based on the configuration interaction approach with a re- 
stricted excitation space like multireference configuration interaction limited to single and 
double excitations (MRCISD) are not size-consistent and therefore they are not well suited 
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for calculations of the interaction energy. 

The most popular example of a one-electron hierarchy is the family of Dunning 
correlation-consistent polarized valence basis sets, cc-pVXZ jl9| with the cardinal num- 
ber X going from D (double-zeta), through T indicating triple- zeta, and so on. These have 



successfully been combined with the HF, MP2, CCSD, CCSD 
models for the calculation of various molecular properties 



20 



'T) hierarchy of wave function 



22|. The basis- set limit, corre- 



sponding to X — )• oo, may be approached eit 
finite cardinal numbers towards infinite X 23 



rer by extrapolating the results obtained with 



24j , or by replacing the standard one-electron 



25 
31 



37fl and shown to accelerate 



38N4Q|. 



hierarchy by explicitly-correlated methods, such as CCSD-F12 and CCSD(T)-F12 
in which the interelectron distance r±2 is explicitly introduced into the wave function 
The F12 methods have recently been implemented efficiently 
the convergence towards the basis-set limit for a number of properties 

In the present paper, we combine all-electron spin-unrestricted CCSD(T)-F12 calcula- 
tions with frozen-core FCI calculations to yield a highly accurate best estimate of the Li-LiH 
interaction potential. We also compare the F12 interaction energies with results obtained 
from standard (not explicitly correlated) CCSD(T) calculations. We then characterize the 
ground-state potential, analyze possible interactions with excited states, and investigate 
channels for reactive collisions. 



II. COMPUTATIONAL DETAILS 

We have calculated the interaction energies between the lithium atom and the lithium 
hydride molecule in Jacobi coordinates (R,r,6), defined for the isotopic combination 7 Li- 
Li 1 H. Calculations were performed for states of 2 A' symmetry in the C s point group. The 
LiH bond distance, r, was initially kept frozen at the LiH monomer equilibrium distance of 



3.014 bohr |4l|. The distance R between Li and the center of mass of LiH ranged from 3.0 
to 10.0 bohr with an interval of 0.5 bohr, and then from 11.0 to 20.0 bohr with an interval 
of 1.0 bohr. Additional distances of 30.0, 40.0, and 50.0 bohr were also used. The angle 
6, between the vector pointing from Li to H in the LiH molecule and the vector pointing 
from the center of mass of the molecule to the Li atom, was varied from 0° to 180° with an 
interval of 15°; 6 = 0° corresponds to Li-H — Li configurations. We thus used a total of 28 
intermonomer distances, R, which combined with the 13 values of 9 yielded 364 grid points 
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on the two-dimensional interaction energy surface. 

Calculations with uncorrelated basis functions were carried out using the unrestricted ver- 
sion of the coupled-cluster model CCSD(T) with Dunning's cc-pVXZ(-mid) basis sets with 
X = D, T, Q, 5, where mid indicates the inclusion of an additional set of basis functions, the 



so-called midbond-95 set 42|, placed at the middle of the Li-LiH distance R. All electrons 
were correlated in these calculations. Additionally, for the purpose of comparison with the 
FCI results (see below), the frozen-core approximation (Ict^h and Isli orbitals kept frozen) 
was used for the cc-pVQZ basis. All these calculations were carried out using the molpro 
package [43] . The full basis set of the dimer was used in the supermolecular calculations and 
the Boys and Bernardi scheme [44j was used to correct for basis-set superposition error 



34 
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The explicitly correlated spin-unrestricted CCSD-F12 and CCSD(T)-F12 
calculations were carried out with the MOLPRO code [43] to establish the CCSD and 
CCSD(T) basis-set limits for the LiH-Li interaction. We chose to use the F12b variant 
35l . |46| of the explicitly correlated spin-unrestricted energy implemented in the MOLPRO 
code. Employing the fixed-amplitude ansatz for the F12 wave function ensured the orbital 
invariance and size-consistency of the CCSD-F12 and CCSD(T)-F12 results. The QZVPP 
basis set j^] was em ploy ed as the orbital basis in the F12 calculations. The corresponding 



QZVPP-ik basis set 

nn 



48] was used as the auxiliary basis for the density-fitting approxima- 
tion 45|, |49| for many-electron integrals, while the uncontracted version of the QZVPP-ik 



basis was used to approximate the Resolution-of-Identity in the F12 integrals |50l . l5l| . In 
addition, the valence correlation in the dimer was described with the full configuration in- 
teraction method (FCI). The FCI and standard CCSD(T) calculations in the frozen-core 

n 

approximation were carried out using the cc-pVQZ basis. The DALTON package [52] and the 
LUCIA program [53] were combined to yield the FCI results. 

To calculate potential energy surface V(R, 9) with the LiH bond length kept fixed at its 
equilibrium value we used computational scheme which was previously applied in theoretical 



54j-|58||. The potential V(R, 9) 



studies of the ground and excited states of the calcium dimer 
was constructed according to the following expression: 

V(R, 9) = v CCSD ^~ F12 (R, 9) + 8V™(R, 9) (1) 

where v CCSD ^- F12 (R, 9) contribution was obtained from all-electron CCSD(T)-F12 calcula- 
tions, while the correction for the valence- valance correlation beyond the CCSD(T)-F12 level, 
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8V™{R, 0), was calculated in an orbital cc-pVQZ basis set. Both terms, \/ CCSD ( T )- F12 ( J R, Q) 
and 8VjSy{R, 0), were obtained from the standard expressions for the supermolecule inter- 
action energy, as given in Ref. js7 |. 

The long-range asymptotic form of the potentials is of primary importance for cold col- 
lisions. We have therefore computed the leading long-range coefficients that describe the 
induction and dispersion interactions up to and including R~ 10 and I = 4 terms, 

10 n-4 

V(R,6) = -Y,J2^( cose ^ ( 2 ) 

n=6 1=0 

where I is even/odd for n even/odd, and C l n = C l n (md)+C l n (disp). The long-range coefficients 
C n (ind) and C^(disp) are given by the standard expressions (see, e.g., Refs. 



59 



). The 



multipole moments and polarizabilities of LiH were computed with the recently introduced 
explicitly connected representation of the expectation value and polarization propagator 



within the coupled-cluster method 



6jJ-|63j, while the Li polarizabilities (both static and 



at imaginary frequencies) were taken from highly accurate relativistic calculations from 
Derevianko and coworkers [64j. 

The interaction potentials were interpolated between calculated points using the repro- 

n 

ducing kernel Hilbert space method (RKHS) [65fl with the asymptotics fixed using the ab 
initio long-range Van der Waals coefficients. The switching function of Ref. \q6\ was used to 
join the RKHS interpolation smoothly with the Van der Waals part in the interval between 
R a = 18 and R b = 26 bohr. 



III. CONVERGENCE OF THE LI LIH INTERACTION POTENTIAL TOWARDS 
THE EXACT SOLUTION 

In sec. IIII Al we analyze the convergence of the Li-LiH interaction potential with respect 
to the one-electron and iV-electron hierarchies. Based on the analysis, we give in sec. IIIIBI 
our best estimate for the ground-state interaction potential with the Li-H bond length fixed 
at its monomer equilibrium value. The features of the potential are presented in sec. IIII CI 
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A. Convergence of the one-electron and iV-electron hierarchies 



In order to investigate the saturation of the Li-LiH interaction energy in the one-electron 
space, we have analyzed three characteristic points of the Li-LiH potential (the global mini- 
mum, the saddle point, the local minimum, and one point very close to the avoided crossing: 
R = 5.5 bohr and 9 = 0.0°). The characteristic points were obtained from the potentials 
calculated at the CCSD(T) / cc-pVXZ-mid level of theory, for X = D, T, Q, and 5. The in- 
teraction energies were then compared to the corresponding energies of the spin-unrestricted 
CCSD(T)-F12 / QZVPP potential (approximation F12b), which serves as the basis-set limit. 
To evaluate the accuracy of the pure one-electron basis (not explicitly correlated), the rel- 
ative percentage errors, Ap 12 b = (y cc ~ pyxz — V F12h ) / |V F12b | • 100%, were determined for 
each X at every characteristic point. The results are given in Tabled 

We have also evaluated the characteristic points from the extrapolated interaction energy 
surfaces, which were generated as follows: at each grid point, the extrapolated total energies 
for Li, LiH, and Li-LiH were obtained by adding the Hartree-Fock energy calculated with 
cardinal number X to the ex trap olated correlation energy, £vx-i)x> obtained from the two- 
point extrapolation formula 



23 



peon- j-icorr i {n\ 

&(x-i)x ~ &x + ^ - - 1' [ } 

where E^ r _^ and E^" are the correlation energies obtained for two consecutive cardinal 
numbers, (X — 1) and X, respectively. The final extrapolated interaction energy at a single 
grid point is obtained by subtracting the Li and LiH extrapolated total energies from the 
Li-LiH extrapolated total energy. We used the values a = 2 and a = 3, which were 



67 



as the ones most 



recommended by Jeziorska et al. in their helium dimer study 
suited for extrapolating all the components of the interaction energy. The energies of the 
characteristic points obtained in this way were compared with the CCSD(T)-F12 / QZVPP 
results and the corresponding values of A^i2b are included in Table [H 

The relative percentage errors, Api2b ; are plotted in Fig. [T] for both plain (non- 
extrapolated) and extrapolated characteristic points. For the global minimum, the plain 
cc-pVXZ results approach the basis-set limit from above and the convergence is smooth and 
fast: the error is reduced by a factor of 2 to 3 for each increment in X. The extrapolation 
accelerates the convergence: the (X — 1)X extrapolated interaction energies have a quality 
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at least that of the plain cc-pV(X + 1)Z results. Though the extrapolation with a = 2 seems 
to be more efficient than that with a = 3 for the DT and TQ cases, it actually overshoots 
the basis-set limit when the Q and 5 cardinal numbers are used. More importantly, using 
a = 2 leads to irregular behaviour: the Q5 extrapolation results in a lower quality than the 
TQ extrapolation. In contrast, extrapolation with a = 3, though slightly less efficient for 
low cardinal numbers, exhibits highly systematic behaviour and leads to an error as small 
as 0.01% for the Q5 extrapolation. 

Similar behaviour of the extrapolation schemes is observed for the point near the avoided 
crossing. Both extrapolations, with a = 2 and a = 3, converge smoothly towards the basis- 
set limit, but the convergence is not as fast as in the case of the global minimum. In contrast 
to the global minimum, there is no problem here with overshooting the basis-set limit. For 
each pair of cardinal numbers (X — 1)X the extrapolation with a = 2 gives results slightly 
more favourable than using a = 3, with the smallest error of 0.19% for the Q5 extrapolation. 

For the saddle point and local minimum, the convergence of the relative errors is not as 
smooth as for the global minimum: the relative error for X=D is surprisingly small. This is 
obviously accidental and does not reflect particularly high quality of the cc-pVDZ basis set. 
Indeed, when the cc-pVDZ results are employed in Eq. ([3]), the extrapolation worsens the 
accuracy: the errors for the DT extrapolation are much larger than the errors for both the 
X=D and X=T plain results, independent of the value of the a extrapolation parameter. 
Starting from X=T, the plain results smoothly approach the basis-set limit, though the 
convergence is clearly slower than in the case of the global minimum. The extrapolation 
with a = 2 is unsystematic and unpredictable, as in the case of the global minimum, while 
that with a = 3 smoothly approaches the basis-set limit. The errors of the Q5 extrapolation 
with a = 3 are —0.49% for the saddle point and —0.13% for the local minimum. 



Patkowski and Szalewicz 



69] recently investigated Ar 2 with the CCSD(T)-F12 method. 



They found that the F12a and F12b variants [35| gave significantly different results. They 
also concluded that, for Ar 2 , calculations with explicitly correlated functions cannot yet 
compete with calculations employing extrapolation based on conventional orbital basis sets. 
Indeed, while their orbital results converged smoothly towards the extrapolated results, the 
CCSD(T)-F12a and CCSD(T)-F12b results behaved erratically with respect to both the 
orbital and the extrapolated results. Table [I] shows that this is not the case for the Li-LiH 
system. In our case the CCSD(T)-F12a and CCSD(T)-F12b results are quite similar and 
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are fully consistent with the plain and extrapolated results with conventional basis sets. It 
should be stressed, however, that Ar 2 is bound mostly by dispersion forces, while the main 
source of the bonding in Li-LiH is the induction energy, which is less sensitive to the basis- 
set quality. This may at least partly explain the success of the CCSD(T)-F12 calculations 
for Li-LiH. 

Finally, it is important to note here that, while the interaction energy at the characteristic 
points varies considerably with the basis set and extrapolation method, the positions of the 
points (i.e., the distance R and angle 9 at which the characteristic points occur) remain 
practically unaffected by the choice of the basis set and extrapolation scheme. 

To analyze the convergence of the CCSD and CCSD(T) models in the iV-electron space, 
Fig. [2] compares the characteristic points (global minimum, saddle point, local minimum, 
and near the avoided crossing) of the Li-LiH potential calculated at the CCSD / cc-pVQZ 
and CCSD(T) / cc-pVQZ levels of theory with the characteristic points obtained at the 
FCI / cc-pVQZ level. The lo"LiH and IsLi orbitals were kept frozen in the calculations. As 
expected, the iV-electron error is reduced by a factor of 3 to 4 when the approximate triples 
correction is included in the calculations. It can also be seen from the figure that the global 
minimum is the most sensitive and the local minimum is the least sensitive to the description 
of the electron correlation. 

B. The best estimate of the ground-state Li LiH potential energy surface 

Because of the negligible one-electron error in the CCSD(T)-F12 calculations and to 
the rather large basis set used in the FCI / cc-pVQZ calculations, and assuming that the 
one-electron and iV-electron errors are approximately independent, the best estimate of the 
ground-state interaction energy surface for the LiH-Li is 

^best = ^CCSD(T)-F12 + Jy FCI + fiyFCl ^ ^ 

where 1/CCSD(t)-fi 2 is the CCSD(T) basis-set limit energy (i.e., the CCSD(T)-F12 result) 
and the FCI correction, 5Vj_*, is obtained by subtracting the CCSD(T) / cc-pVQZ energy 
from the FCI / cc-pVQZ energy, both calculated in the frozen-core approximation. The 
quantity 5V FCl accounts for the last remaining correction (in the non-relativistic limit), 
namely the effects of core-core and core-valence correlation in the FCI / cc-pVQZ calcula- 
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tions, 



5V™ = 5V™ n - 6V™, (5) 

where the subscript "all" refers to all electrons correlated. 

The quantity 5V FC1 is a measure of the uncertainty in our best estimate V hest . To estimate 
this, we may safely assume that SV FG1 is at most as large as the corresponding 5V^ T \ 

5V™ < SV™ = 5V^L - 5V™, (6) 

where 

JfT/(T) _ T/ CCSD(T) T/C CSD f 7 \ 

V all-all — "all-all "all-all I ' I 

*v£3 = K C -? D(T) - e c v SD , (8) 

with V^J^jj , Kn-al?' y^-v D ^ T \ an d K-v SD denoting interaction energies calculated at 
the CCSD(T) / cc-pVQZ or CCSD / cc-pVQZ level, correlating all electrons or using the 
frozen-core approximation, as appropriate. As can be seen from Fig. [21 the differences 
between CCSD(T) and CCSD are, for the characteristic points of the potential, 2 to 3 times 
larger (and for the rest of the potential at least 1.5 times larger) than the differences between 
FCI and CCSD(T). Eq. ([HD is therefore actually a conservative estimate for 5V FCl . The root 
mean square error for 5V^ /5V^}_ all , over the whole potential is 4.1%. We thus consider 
that our best estimate of the ground-state interaction energy for LiH-Li, Eq. (j3j), has a 
(conservative) total uncertainty of 5% of the FCI correction (5Vj_S^). The analysis of the 
Li-LiH potential in the remainder of this paper is based on the interaction energies obtained 
using Eq. unless otherwise stated. 

To justify our error estimation we have performed calculations with all electron correlated 
at the FCI level for the set of characteristic points of the potential. Due to the immense 
memory requirements of the FCI calculations with seven electrons we were able to apply the 
cc-pVDZ basis set only. The FCI/cc-pVDZ results together with the CCSD(T)/cc-pVDZ, 
both with and without the frozen-core approximation, are presented in Table HT1 The error in 
the FCI correction calculated with frozen core is as small as 0.76 % for the examined points. 
We may see that the approximation with the FCI valence correction added to CCSD(T), Eq. 
([I)), reproduces the exact FCI results with accuracy better than 1% of the FCI correction 
(SVjS^). This confirms our estimate of 5% uncertainty in the FCI correction 8V F _y- 
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C. Features of the ground-state potential energy surface 



In Table HTT1 we have listed the characteristic points of the potential energy surfaces of the 
ground state, which correlates at long range with Li( 2 S) + LiH (X 1 E+), and the first excited 
state, which correlates at the long range with Li( 2 P) + LiH (X 1 E + ). Both these states are 
of 2 A' symmetry in the C s point group. The latter is included in Table II III since, as will be 
discussed in the next section, it shows an avoided crossing with the ground-state potential 
for the linear LiH-Li geometry. Table [LTTl shows that the interaction potential for the ground 
state of Li-LiH is deeply bound, with a binding energy of 8743 cm -1 at the global minimum. 
The global minimum is located at a skew geometry with i? e =4.40 bohr and # e = 46.5°, and 
is separated by a barrier around _R=6.3 bohr and 0=136.0° from a shallow local minimum 
at the linear Li-LiH geometry. The local minimum is at _R=6.56 bohr, with a well depth 
of only 1623 cm -1 . The excited-state potential shows only one minimum, at _R=5.66 bohr, 
with a binding energy of 4743 cm -1 . 

A contour plot of the ground-state potential is shown in the left-hand panel of Fig. [31 while 
the full-CI correction to the CCSD(T) potential, 5Vj^, is shown in the right-hand panel. 
The correction is very small compared to the best potential. It amounts to 0.4% around 
the global minimum, and approximately 1% at the local minimum. Thus, our estimated 
error of the calculation, 5% of the full-CI correction, translates into 0.05% error in the 
potential itself. We would like to reiterate here that such a small error was achieved not 
only because the interelectron distance was included explicitly in the ab initio CCSD(T)- 
F12 calculations, but also because of the very small valence- valence correlation beyond the 
CCSD(T) level. The smallness of the valence-valence correlation beyond the CCSD(T) level 
is not so surprising, since Li-LiH has only three valence electrons, and the exact model for 
a three-electron system would be CCSDT, coupled-cluster with single, double, and exact 



triple excitations 



701 ] . Our results show that the triples contribution to the correlation 



energy beyond the CCSD(T) model for the valence electrons is very small. 

The potential for the ground state of Li-LiH is very strongly anisotropic. This is easily 
seen in the left-hand panel of Fig. |3l and in Fig. HI which shows the expansion coefficients 
of the potential in terms of Legendre polynomials Pj(cos0), 

oo 

V(R,9) = ^2v l (R)P l (cos9). (9) 

1=0 
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Here, Vq(R) is the isotropic part of the potential and {V;(i2)}^ 1 is the set of anisotropic 
coefficients. Fig. H] shows that, around the radial position of the global minimum, i?=4.36 
bohr, the first anisotropic contribution to the potential, Vi(R), is far larger than the isotropic 
term, Vq(R). The higher anisotropic components, with / = 2,3, etc., contribute much less 
to the potential. 

As mentioned above, calculations of collision dynamics at ultralow temperatures require 
accurate values of the long-range potential coefficients, Eq. (j2J). Some important scattering 
properties, such as the mean scattering length and the heights of centrifugal barriers, are 
determined purely by the Van der Waals coefficients. The calculated coefficients for Li- 
LiH are presented in Table IIVI Because of the large dipole moment of lithium hydride 
and the relatively high polarizability of the lithium atom, the lowest-order, most important, 
coefficients are dominated by the induction contribution. For example, the induction part 
of Cl and C\ is 887 a.u., which accounts for 71% of C° and 98% of C\. 

IV. INTERACTION BETWEEN THE GROUND AND EXCITED STATES 

A. Low-lying excited state potential, nonadiabatic coupling matrix elements, and 
diabatic potentials 

We encountered convergence problems with CCSD(T) calculations at the linear LiH- 
Li geometry around R=5.6 bohr, due to the presence of a low-lying excited state. The 
excited state correlates with the Li( 2 P)+LiH(X 1 S + ) dissociation limit, but closer investiga- 
tion revealed that, at linear Li-HLi geometries near the crossing with the ground state, it 
has ion-pair character, Li + ( 1 S) + LiH~( 2 £). The ion-pair state itself has a crossing near 
R = 9 bohr with the lowest 2 A' state correlating with Li( 2 P) + LiH(X 1 S). This is shown 
schematically in Fig. [5j Away from linear Li-HLi geometries, the excited state has covalent 
character and remains below the ion-pair state all the way to dissociation. The avoided 
crossing between the ground state and the first excited state is at i?=5.66 bohr, which is 
near the minimum of the ground-state potential at the linear geometry, and the energetic 
distance between the two states at the avoided crossing is only 94 cm -1 . 

In order to investigate how far the excited state may affect the scattering dynamics, 
we computed the full potential energy surface for the excited state in question by means 
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of equation-of-motion coupled-cluster method with single and double excitations (EOM- 
CCSD) 7lM73| implemented in the QCHEM code [74], using the orbital cc-pVQZ basis set. 
Cuts through the ground-state and excited-state potential energy surfaces at selected values 
of the angle 9 are shown in Fig. |6j It may be seen that it is only near the linear LiH- 
Li geometry that the two states come very close together. If we distort the system from 
the linear geometry, the excited state goes up in energy very rapidly, and around the global 
minimum energy, 9 « 45°, it is almost 6000 cm -1 above the ground state. The importance of 
the possible interaction between the ground and excited states can be measured by analyzing 
the (vectorial) nonadiabatic coupling matrix elements T12, defined as T12 = (^iIV^)) where 
V is the gradient operator of the position vector and ^1 and ^2 are the wave functions 
of the two lowest states. On the two-dimensional surface, we may define radial t\2,r = 
(^i\d^2/dR) and angular Ty2.fi = (^i\d^2/d9) components of the vector T\2- We evaluated 
T12 for all (R, 9) geometries by means of the multireference configuration interaction method 



limited to single and double excitations (MRCI) [75|, l76| . using the molpro code |43j. The 
nonadiabatic coupling is largest when the two states are very close in energy, as it can be seen 
in Fig. [6j While at 9 = 0° the radial component of the nonadiabatic coupling approaches 
the Dirac delta form near the crossing point i? ac , with increasing angle it becomes a broad 
function of approximately Lorentzian shape. [Note the different scales on the vertical axes 
of the different panels.] 

The transformation from the adiabatic representation to a diabatic representation may 
be expressed in terms of a mixing angle 7, 

Hi = V 2 sin 2 7 + V\ cos 2 7, H 2 = V\ sin 2 7 + V 2 cos 2 7, H 12 — (V2 — Vi) sin7C0S7, (10) 

where V\ and V2 are the ground-state and excited-state adiabatic potentials, Hi and H2 are 
the diabatic potentials, and if 12 is the diabatic coupling potential. In principle, the mixing 
angle 7 may be obtained by performing line integration of the nonadiabatic coupling T\%, 

7(^=7(^0)+ [ R Ti 2 dl, (11) 

where Ro is the starting point of the integration. For polyatomic molecules, however, the 
mixing angle 7 obtained by integrating this equation is non-unique due to the contributions 
from higher states. To circumvent the problem of path dependence, one may assume that 
we deal with an ideal two-state model. 
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In our case, however, the ion-pair surface Li + ( 1 S)+LiH~( 2 E) shows another crossing at 
small angles and large distances, 9 < 15° and R ~ 9 bohr, with another excited-state 
potential that correlates with the Li( 2 P)+LiH( 1 S) dissociation limit. Thus a third state \l/3 
comes into play and a two-state model is not strictly valid. The energy of the first excited 
state goes up very rapidly with the angle 9, and at the same time the contribution of the 
ion-pair configuration to the wave function of the first excited state, diminishes rapidly. 
Fortunately, the nonadiabatic coupling matrix elements between the two lower states r 12 
and between the two higher states T23 are well isolated. The maximum of r 12 is separated 
from the maximum of T23 by more than 4 bohr; the locations of the crossing points between 
the surfaces for 9 = are shown in Fig. [5j Moreover, the coupling T13 between the ground 
state and the third state is negligible over the whole configurational space. Thus, following 



the discussion of Baer et al. 



77| on the application of the two-state model, we conclude that 



the necessary conditions are fulfilled for the Li-LiH system. Due to the spatial separation 
of the nonadiabatic couplings T12 and T23, using the diabatization procedure based on the 
two-state model is justified. It is worth noting that in our particular case we could not use 
the so-called quasi-diabatization procedure [78], since it is not possible to assign a single- 
reference wave function. This is due to the fact that the excited state shows admixture from 
the ion-pair state. 

As the starting point of the integration in Eq. fill I) , we chose R = 20 bohr and 9 = 0° 
and followed a radial path along 9 = 0° and subsequently angular paths at constant R. 
The diabatic potentials were then generated according to Eqs. (JTU]). Contour plots of the 
adiabatic, diabatic, and coupling potentials, and of the mixing angle 7, are presented as 
functions of R and 9 in Fig. [7J We consider first the mixing angle 7, which is plotted in the 
bottom right-hand panel of Fig. [7J As expected, the mixing angle shows an accumulation 
point at 9 = 0° at a distance R corresponding to the closely avoided crossing between the 
ground and excited states. For 9 = 180°, the mixing angle is non-negligible, even at large 
distances. The coupling potential if 12 vanishes quite slowly with distance R, as R~ 3 . For 
the coupling between the ground and ion-pair states, this long-range decay is exponential, 
because of the different dissociation limits of the two surfaces. As expected, at large distances 
the two diabatic surfaces approach the respective adiabatic surfaces. The diabatic surface 
that correlates asymptotically with the excited-state Li( 2 P)+LiH surface has an important 
contribution from the ground-state adiabatic potential only inside the avoided crossing and 
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at small angles 9. The diabatic surface that correlates asymptotically with the ground state 
resembles the ground-state adiabatic surface rather less closely, especially at large values 
of 9. The coupling between the diabatic states is small over a significant region of 9 and 
LiH bond length r in the vicinity of the crossing. Physically, this means that the dynamics 
will be strongly nonadiabatic in this region, and to take this rigorously into account would 
require a full two-state treatment of the dynamics. However, there are no open channels 
that involve the second surface, and any collisions that cross onto it must eventually return 
to the original surface. Its effect in collision calculations will therefore be at most to cause 
a phase change in the outgoing wavefunction. 



B. Conical intersection 



It is well known that potential ener gy s urfaces for homonuclear triatomic systems com- 
posed of hydrogen j^l or lithium atoms (80 1 show conical intersections at equilateral triangu- 
lar geometries. Analogous behaviour may be expected for Li 2 H, at geometries where the two 
lithium atoms are equivalent, i.e., C 2v geometries. Thus far, our discussion of the potential 
for Li-LiH has been restricted to two dimensions with the bond length of the LiH molecule 
fixed at its equilibrium value, and no conical intersection was observed. However, if we start 
to vary the bond length of the LiH molecule, conical intersections show up immediately. 

At C 2v geometries, with the two LiH bond lengths equal, there are two low-lying electronic 
states, of 2 Ai and 2 B 2 symmetries, that cross each other as a function of the internuclear 
coordinates. Fig. [8] shows contour plots of the two potential energy surfaces and of the 
difference between them, and the top panel of Fig. O summarizes some key features of the 
surfaces. The 2 Ai state has a minimum energy of —8825 cm -1 at r(LiH) = 3.22 bohr and 
an Li-H-Li angle of 95°. MRCI calculations with all coordinates free to vary confirm that 
this is indeed the absolute minimum geometry. There is also a saddle point on the 2 Ai 
surface at a linear H-Li-H geometry with r(LiH) = 3.04 bohr and an energy of —4992 cm -1 , 
which is a minimum in symmetry. The 2 B 2 state has a minimum energy of —5136 

cm -1 at r(LiH) = 3.17 bohr at a linear Li-H-Li geometry. The 2 Ai saddle point and 2 B 2 
linear minimum have symmetries 2 E+ and 2 £+ respectively in -Dooh symmetry, but mix and 
distort if the constraint on the LiH bond lengths is relaxed, to form a 2 S + state in Coo V 
symmetry with a minimum at a linear geometry with r(LiH) distances of 3.00 and 3.33 bohr 
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and an energy of —5323 cm -1 . Even this is a saddle point with respect to bending on the 
full potential surface in C s symmetry. 

The 2 Ai and 2 B2 states are of different symmetries at C*2 V geometries, but both are of 
2 A' symmetry when the geometry is distorted from C*2 V to C s symmetry. The two states 
therefore mix and repel one another at geometries where the two LiH bond lengths are 
different, but a seam of conical intersections runs along the line where the energy difference 
is zero at geometries. 

The fixed LiH distance used in previous sections (r = 3.014 bohr, shown as a dashed 
line on the figure) keeps the 2 Ai surface just below the 2 B2 surface. However, if we allow 
for the vibrations of LiH, the seam of conical intersections becomes accessible at near-linear 
LiH-Li geometries, where the zero of the energy difference appears for an Li-H distance only 
slightly larger than 3.014 bohr. At non-linear geometries the seam quickly moves to Li-H 
distances far outside the classical turning points of the ground vibrational level of free LiH, 
which are 2.72 and 3.35 bohr. 

It is interesting to compare the features of the conical intersections in Li2H with those in 
other triatomic molecules formed from Li and H atoms: LiH2, Li3 and H3. In the case of 
LiH 2 , the seam of intersections occurs at highly bent geometries with an angle between 
the two Li-H bonds of approximately 30° and arises from degeneracy between the surfaces 
of Ai and B 2 symmetry. The global minimum of B 2 symmetry is located at r(LiH) = 3.23 
bohr and a bond angle H Li H of 28° |8l|. This contrasts with Li2H, where the minimum 
of B2 symmetry is at a linear Li-H-Li configuration. The energy of the lowest point on the 
seam of intersections is about 9000 cm -1 above the Li( 2 S)+H2(A 1 S 9 ) threshold, so that it 
is irrelevant for low and medium-energy collisions between Ho and Li in their ground states, 



though it is important for quenching of Li( 2 P) by H 2 8l|, |82|. 

The conical intersections for the doublet states of Li 3 and H 3 occur at equilateral trian- 
gular geometries, where the ground state is doubly degenerate and has symmetry 2 E' in the 
Dm point group. In the case of H3, the lowest-energy point on the seam is located at an 
energy more than 20000 cm -1 above the H( 2 S)+H2(A 1 S 9 ) threshold, so that nonadiabatic 



effects are negligible in t 
geometric phase effects 



3 2 collisions 



83] , although the conical intersection also produces 



84]. For Li 3 , the energetics are essentially different. The lowest- 

ow the Li( 2 S)+Li 2 (A 1 S 5 ) threshold and 



)C 



energy point on the seam is around 4000 cm 1 
only 500 cm -1 above the C^ v global minimum 85|. This is likely to produce considerable 
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nonadiabacity in collisions of Li 2 with Li. 

To conclude, in all the triatomic molecules formed from H and Li there are seams of 
crossings that occur at configurations of the highest possible symmetry, either C% v or D^. 
For Li3 and Li2H the conical intersections are accessible during atom-molecule collisions, 
while for H3 and LiH2 nonadiabatic processes are unimportant if the colliding partners are 
in their ground states and have relatively low kinetic energy. 



V. REACTION CHANNELS 

Several reaction channels exist that might affect sympathetic cooling [86] in Li+LiH. 
These are the exchange reaction, 

LiH + Li -> Li + HLi (12) 

and two insertion reactions, 



'1 



LiH + Li -> Li 2 + H, (13) 

producing Li2(X 1 E+) and Li 2 (a 3 £+) plus a ground-state H atom. The energetic location of 
the entrance and exit channels of these reactions, as well as those of the potential minima for 
linear and C 2v geometries, are shown in the upper panel of Fig. |9j The insertion reactions 
are highly endothermic, with an energy difference between the entrance and exit channels of 
the order of 12000 cm" 1 and 22500 cm" 1 for Li 2 (X 1 S+)+H and Li 2 (a 3 £+)+H, respectively. 

To make the discussion more quantitative, Fig. O also shows two-dimensional plots of the 
energy as functions of the internal coordinates. For the exchange reaction, we held Li-H-Li 
at linear geometries and varied the distances from the two lithium atoms to the hydrogen 
atom. For the insertion reaction, Li-Li-H was kept bent, with the angle Z(HLilLi2) held 
constant at the C 2v equilibrium value 42.5°, while the Li-Li and Li H distances were varied. 
To make the plots consistent with the correlation diagram shown on the upper panel, the 
zero of energy was fixed at that of Li-LiH separated to infinite distance with the Li-H bond 
length fixed at the monomer equilibrium value. 

Let us consider the exchange reaction first. The two-dimensional cut through the potential 
energy surface is presented in the left-hand panel of Fig. [9j The potential energy surface 
of linear Li 2 H has two equivalent minima with an energy of —5323 cm -1 , separated by a 
small barrier 187 cm^ 1 high. The linear minima are in any case substantially above the 
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absolute minimum (8825 cm -1 ), so this small barrier will have no important effect on the 
collision dynamics. The exchange reaction produces products that are indistinguishable 
from the reactants, so reactive collisions cannot be distinguished from inelastic collisions 
experimentally (unless the two Li atoms are different isotopes). 

An analogous two-dimensional cut through the potential energy surface corresponding to 
Li 2 (X 1 E 9 )+H products is presented in the right-hand panel of Fig. The plot illustrating 
the reaction to form Li 2 (a 3 E u )+H products is not reported, as the reaction is even more 
endothermic. The surface includes the absolute minimum at an energy of —8825 cm -1 . The 
entrance channel for this reaction corresponds to an Li-H distance of 3.014 bohr at large 
Li-Li distance, while in the exit channel the Li-Li distance is approximately 5.05 bohr when 
the Li-H distance is very large. However, this reaction cannot occur at low collision energies. 

VI. SUMMARY AND CONCLUSIONS 

In the present paper, state-of-the-art ab initio techniques have been applied to compute 
the ground-state potential energy surface for Li-LiH in the Born-Oppenheimer approxima- 
tion. The interaction potential was obtained using a combination of the explicitly-correlated 
unrestricted coupled-cluster method with single, double, and approximate noniterative triple 
excitations [UCCSD(T)-F12] for the core-core and core-valence correlation, with full config- 
uration interaction for the valence-valence correlation. The main results of this paper can 
be summarized as follows: 

1. The Li-LiH system is strongly bound: if the LiH bondlength is held fixed at the 
monomer equilibrium distance of 3.014 bohr, the potential energy surface has a global 
minimum 8743 cm -1 deep at a distance i?=4.40 bohr from the lithium atom to the 
center of mass of LiH, and a Jacobi angle 9 = 46.5°. It also shows a weak local 
minimum 1623 cm -1 deep at the linear Li-LiH geometry for i?=6.56 bohr, separated 
from the global minimum by a barrier at i?=6.28 bohr and 9 = 136°. If the LiH bond 
length is allowed to vary, the potential minimum is at a depth of 8825 cm -1 , at a C 2v 
geometry with LiH bond length of 3.22 bohr and an Li-H-Li angle of 95°. 

2. The full-CI correction for the valence-valence correlation to the explicitly correlated 
CCSD(T)-F12 potential is very small. The remaining error in our calculations is due 
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to the neglect of the core-core and core-valence contributions, and is estimated to be 
of the order of 0.05% of the total potential. 

3. To evaluate the performance of the conventional orbital electron-correlated methods, 
CCSD and CCSD(T), calculations were carried out using correlation-consistent polar- 
ized valence X-tuple zeta basis sets, with X ranging from D to 5, and a very large 
set of mid-bond functions. Simple two-point extrapolations based on the single-power 
laws X~ 2 and X~ 3 for the basis-set truncation error reproduce the CCSD(T)-F12 
results for the characteristic points of the potential with an error of 0.49% at worst. 

4. The potential for the ground state of Li-LiH is strongly anisotropic. Around the 
distance of the global minimum, the isotropic potential Vo(R) is almost two times 
smaller than the first anisotropic contribution V\ (R). Higher anisotropic components, 
with / = 2, 3, etc., do not contribute much to the potential. 

5. At the linear LiH-Li geometry, the ground-state potential shows a close avoided cross- 
ing with the first excited-state potential, which has ion-pair character around the 
avoided crossing point. The full potential energy surface for the excited state was ob- 
tained with the equation-of-motion method within the framework of coupled-cluster 
theory with single and double excitations. The excited-state potential has a single 
minimum 4743 cm -1 deep for the linear LiH-Li geometry at i?=5.66 bohr. The en- 
ergy difference between the ground and excited states at the avoided crossing is only 
94 cm -1 . An analysis of the nonadiabatic coupling matrix elements suggests that 
dynamics in the vicinity of the avoided crossing will have nonadiabatic character. 

6. When stretching the LiH bond in the Li-LiH system, a seam of conical intersections 
appears for C*2 V geometries, between the ground state of 2 Ai symmetry and an excited 
state of 2 B2 symmetry. At the linear LiH-Li geometry, the conical intersection occurs 
for an Li-H distance which is only slightly larger than the equilibrium distance of the 
LiH monomer, but for significantly non-linear geometries it moves to Li-H distances 
far outside the classical turning points of LiH. 

7. The Li-LiH system has several possible reaction channels: an exchange reaction to 
form products identical to the reactants, and two insertion reactions that produce 
Li 2 (a 3 S^) and Li 2 (X 1 £+) plus a ground-state hydrogen atom. The insertion reactions 
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are highly endothermic, with the energy difference between the entrance and exit chan- 
nels of the order of 12000 cm" 1 and 22500 cm" 1 for Li 2 (X 1 S+)+H and Li 2 (a 3 £+)+H, 
respectively. 



In a subsequent paper 87| we will analyze the dynamics of Li-LiH collisions at ultralow 
temperatures, based on our best ab initio potential. We will analyze the impact of the present 
inaccuracies in the ab initio electronic structure calculations, and discuss the prospects of 
sympathetic cooling of lithium hydride by collisions with ultracold lithium atoms. 
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TABLE I: Performance of various orbital basis sets and extrapolation schemes compared to the 
CCSD(T)-F12a and CCSD(T)-F12b results at the characteristic points of the ground state potential 
energy surface of Li-LiH. The notation VXZ with X=D, T, Q, and 5 denotes the result of the orbital 
CCSD(T) calculations in the cc-pVXZ basis with midbond-95 set, XYa with X and Y = D, T, Q, 
and 5, and a = 2 or 3 denote the extrapolated result according to Eq. ([3]), while F12a and F12b 
stand for the explicitly correlated CCSD(T) results with the a and b approximation schemes. Api2a 
and Api2b are the percent error of given result with respect to the CCSD(T)-F12a and CCSD(T)- 
F12b results, respectively. 

V (cm -1 ) A F i2a A F12b V (cm -1 ) A F i 2a A F i 2b V (cm -1 ) A F i 2a A F i 2b V (cm -1 ) A Fi2a A Fi2b 



Global Minimum Saddle point Local minimum Avoided crossing 



VDZ -8547.77 1.87 1.84 
VTZ -8652.38 0.67 0.64 
VQZ -8683.85 0.31 0.27 
V5Z -8698.84 0.14 0.10 


-1548.64 -0.39 -0.61 
-1566.01 -1.52 -1.73 
-1555.93 -0.86 -1.08 
-1551.91 -0.60 -0.82 


-1616.54 -0.65 -0.77 
-1629.62 -1.46 -1.59 
-1618.22 -0.75 -0.88 
-1612.70 -0.41 -0.53 


-4771.50 4.59 4.65 
-4892.57 2.17 2.23 
-4955.39 0.91 0.97 
-4974.35 0.53 0.60 


DT2 -8683.31 0.32 0.28 
TQ2 -8704.14 0.08 0.04 
Q52 -8714.77 -0.05 -0.08 


-1575.59 -2.14 -2.36 
-1541.85 0.05 -0.17 
-1543.21 -0.04 -0.25 


-1637.22 -1.94 -2.06 
-1601.47 0.29 0.17 
-1601.76 0.27 0.15 


-4934.30 1.33 1.40 
-4983.15 0.36 0.42 
-4991.70 0.19 0.25 


DT3 -8668.66 0.48 0.45 
TQ3 -8695.37 0.18 0.14 
Q53 -8707.24 0.03 0.01 


-1571.05 -1.84 -2.06 
-1547.94 -0.35 -0.56 
-1546.78 -0.27 -0.49 


-1633.62 -1.71 -1.84 
-1608.71 -0.16 -0.28 
-1606.24 -0.01 -0.13 


-4914.53 1.73 1.79 
-4971.14 0.60 0.66 
-4984.59 0.33 0.39 


F12a -8710.85 0.00 -0.04 
F12b -8707.77 0.04 0.00 


-1542.60 0.00 -0.21 
-1539.31 0.21 0.00 


-1606.14 0.00 -0.12 
-1604.17 0.12 0.00 


-5001.10 0.00 0.06 
-5004.15 -0.06 0.00 
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TABLE II: Performance of the valence-valence FCI correction SV^Py against exact FCI results 
for the characteristic points of the LiH-Li potential (in cm -1 ). Calculations were done in the 
cc-pVDZ basis set. Subscript all-all refers to all electrons correlated, while v-v denotes frozen- 
core results. A is the percentage error of the SV^S} approximation with respect to the SV^}^: 
Sv FCl -8v FC1 





GM 


SP 


LM 


AC 


^CCSD(T) 


-7553.16 


-1407.50 


-1502.25 


-3785.91 


yFCl 

' v— V 


-7587.86 


-1437.27 


-1522.05 


-3806.93 


5V FCl 


-34.70 


-29.77 


-19.80 


-21.02 


T/ CCSD(T) 
•all-all 


-7590.40 


-1415.56 


-1509.04 


-3842.01 


yFCl 
"all-all 


-7625.07 


-1445.32 


-1528.78 


-3863.46 


"all-all 


-34.67 


-29.76 


-19.74 


-21.45 


Af%l 


-0.09 


-0.04 


-0.30 


0.76 
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TABLE III: Characteristic points of the interaction potentials for the ground state Li( 2 S) + LiH 
(X 1 E + ) and the first excited state, which correlates asymptotically with Li( 2 P) + LiH (X 1 S + ). 



R [bohr] 6 [degrees] V [cm -1 ] 

Ground state 
Global minimum 4.40 
Local minimum 6.56 
Saddle point 6.28 

Excited state 
Global minimum 5.66 0.0° -4743 



46.5° -8743 
180.0° -1623 
136.0° -1565 
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TABLE IV: Long-range coefficients (in atomic units) for Li-LiH, from perturbation theory. Num- 
bers in parentheses indicate powers of 10. 



I -> 1 2 3 4 

C\ 1247.8 869.7 

C\ 8304.2 2902.3 

C\ 4.83(4) 4.24(4) 8923.0 

C\ 3.94(5) 6.42(4) 

C[ 2.03(6) 2.02(6) 2.44(5) 
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Global minimum Saddle Point Local Minimum 





FIG. 1: Relative percentage errors of the interaction energy at the characteristic points (global 
minimum, saddle point, local minimum, and near the avoided crossing: R—5.5 bohr and 6 = 0°) 
of the LiH-Li potential calculated at the CCSD(T) / cc-pVXZ-mid level, where X = D, T, Q, 
5 and mid stands for the midbond-95 set. The errors for the characteristic points obtained by 
extrapolating the plain basis-set results with the two-point extrapolation formula are also shown 
for a = 2 and a = 3. The errors were obtained by comparison with the CCSD(T)-F12b / QZVPP 
results. 
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FIG. 2: The TV-electron error of the characteristic points of the LiH-Li interaction potential calcu- 
lated at the CCSD / cc-pVQZ and CCSD(T) / cc-pVQZ levels of theory. The error was determined 
by comparison with the FCI / cc-pVQZ interaction potential. The l<TLiH l^Li frozen-core approxi- 
mation was used. 
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FIG. 3: Contour plots of the best ab initio potential for the ground state of Li-LiH (left-hand 
panel), and of the full-CI correction to the CCSD(T)-F12 potential (right-hand panel). Energies 
are in cm -1 . 



32 



6000 



4000 



2000 







3 -2000 - 



-4000 - 



-6000 




8 



3 4 5 6 7 

R [bohr] 

FIG. 4: The Legendre components Vi(R) {I = 0, 1, 2, 3, 4) of the ground-state Li( 2 S) + LiH (X X S+) 
interaction potential, see Eq. 
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FIG. 5: Ground-state and excited-state potentials at the linear geometry LiH-Li: demonstration 
of the ion-pair nature of the excited-state potential. 
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FIG. 6: Cuts through the potential energy surfaces for the ground and the first excited states of 
2 A' symmetry for selected values of the angle 6, and the corresponding radial nonadiabatic coupling 
matrix elements as functions of the distance i??5 
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FIG. 7: Adiabatic potentials (top two panels), diabatic potentials (middle two panels), and coupling 
potential (bottom left-hand panel), and the mixing angle 7 (bottom right-hand panel), as functions 
of the geometrical parameters R and 9. Energies are in cm -1 . 
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FIG. 8: Potential energy surfaces for the 2 Ai (top left-hand panel) and 2 E>2 (bottom left-hand 
panel) states of L^H in Civ symmetry, and the difference between them (top right-hand panel), 
which is zero along the seam of conical intersections. Also shown is the coordinate system used for 
C2v geometries. Energies are in cm -1 . 
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FIG. 9: Schematic representation of the possible reaction pathways in collisions of the lithium atom 
with the lithium hydride molecule (upper panel), and two-dimensional cuts of the reactive poten- 
tial energy surfaces for the exchange (left-hand panel) and insertion (right-hand panel) reactions. 
Energies are in cm -1 . 
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